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I. INTRODUCTION 


Although the LHC has unearthed the standard model (SM) Higgs boson |IH3J, so far, it has 
failed to identify any sign of new physics beyond the SM. New physics at the TeV scale is invoked 
in the form of supersymmetry, extra dimensions etc. to address the issue of the hierarchy problem. 
Absence of such new physics can imply that the solution of the hierarchy problem lies somewhere 
else. Hence, it is likely that the SM could be valid up to the Planck scale (Mpi), above which 
gravity is expected to dominate over other forces. 

However, we know that the presence of dark matter (DM) and nonzero neutrino masses indicate 
physics beyond the SM. Hence, it is worthwhile to explore a scenario where the SM, augmented 
by a low energy DM model, is valid up to M P \. Earlier we had explored such a scenario where a 
gauge singlet scalar is added to the SM to pose as a weakly interacting massive particle (WIMP) 
DM [4j. In the present work, we extend the SM by a scalar doublet protected by a Z 2 symmetry. 
The model is popularly known in the literature as the inert doublet (ID) model, first proposed by 
Deshpande and Ma [5]. 

Given the experimental measured values for M t , Mh, and a s , the clectroweak (EW) vacuum is 
reported to reside in a metastable state within the framework of the SM. This has been verified in 
the literature [shui] only very recently using state of the art next-to-next-to leading order (NNLO) 
loop corrections contributing to the Higgs effective potential. In Ref. [T], this analysis was extended 
to the case of a singlet scalar DM model. The stability of the EW vacuum was shown to depend 
on new physics parameters. In this paper, we extend such an analysis to the ID model. We assume 
that ID DM is the only DM particle which saturates the entire DM relic density. In this context, 
we review the constraints on the parameters of the ID model. 

A detailed study on the ID parameter space was recently performed in Refs. [min] indicating 
bounds from EW stability, perturbativity, collider study, clectroweak precision tests (EWPT), 
etc., considering DM annihilation to two-body final states only. In Refs. [HI [13] DM relic density 
was calculated including three-body final states as well. An updated detailed analysis on the 
ID parameter space was performed in Refs. mm reflecting the impact of the M ^ & 126 GeV 
Higgs boson discovery at the LHC. The authors of Refs. mm looked at renormalization group 
evolution (RGE) of model parameters to check the validity of the ID model at higher energies. 
The constraints from the measurements of the diphoton decay channel of the Higgs boson on the 
ID parameter space was also discussed. High scale validity of this model in the presence of a 
right-handed neutrino has been looked at in Ref. I16j . The influence of ID in the Higgs effective 
potential was studied in Ref. ra to explore the possibility of electroweak symmetry breaking a la 
Coleman-Weinberg ra. High scale validity of the two Higgs doublet model has been explored in 
Refs. P2J, :2D] with broken Z 2 , symmetry considering tree level potential only. 

In this paper we improve earlier studies on ID model parameter space by including radiative 
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corrections to the scalar potential to explore the stability of the electroweak vacuum. In particular, 
we permit the Higgs quartic coupling to assume small negative values to render a metastable EW 
vacuum. As our analysis improves the scalar potential considering radiative corrections, the tree 
level stability constraints used in earlier analyses need to be reviewed. In addition, as we demand 
that the theory be valid up to Mpi, the model parameters at the EW scale are constrained by 
the requirement that they satisfy all the bounds up to the M-p\. If in the near future the direct 
detection experiments and colliders confirm the ID model as the dark matter model realized in 
nature, then our study will help to estimate the lifetime of the electroweak vacuum, especially if 
it still remains in the metastable state. As metastability of EW vacuum is considered to be an 
indication of new physics, if the vacuum is metastable even after consideration of the ID model, 
then that could indicate the presence of additional new physics before Mp\. 

The paper is organized as follows. Section |TT| starts with an introduction of the ID model, 
followed by a discussion of the effective scalar potential and renormalization group (RG) running of 
parameters of the model. Constraints on the model are discussed in Sec. HI Tunneling probability 
of the metastable vacuum and the restrictions on the parameters to avoid a potential unbounded 
from below, are mentioned in Sec. IV A detailed study of the parameter space identifying regions 
of EW vacuum stability and metastability is performed in Sec. [V] with the help of various phase 
diagrams. Section VI contains a short discussion on Veltman’s conditions in the context of the ID 


model. We finally conclude in Sec. VII 


II. INERT DOUBLET MODEL 


In this model, the standard model is extended by adding an extra SU( 2) doublet scalar, odd 
under an additional discrete Z 2 symmetry. Under this symmetry, all standard model fields are 
even. The Z 2 symmetry prohibits the inert doublet to acquire a vacuum expectation value. 

The scalar potential at the tree level is given by 


U(d>i, <f> 2 ) = ^|$i| 2 + A^il 4 + p||d> 2 | 2 + A 2 |<h 2 | 4 

+A 3 |$i| j |$ 2 | 2 + A 4 | < hJ < I ) 2 |“ H—— ( < hj$ 2 )' + H.c. 
where the SM Higgs doublet dq and the inert doublet <f> 2 are given by 


( 2 . 1 ) 


( 72 ( V + k + iG ° ) ) ’ 2 ( 75 ( H + ) 

G ± and G° are Goldstone bosons and h is the SM Higgs. 

<f> 2 contains a CP even neutral scalar H , a CP odd neutral scalar A, and a pair of charged 
scalar fields H ± . The Z 2 symmetry prohibits these particles to decay entirely to SM particles. The 
lightest of H and A can then serve as a DM candidate. 
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After electroweak symmetry breaking, the scalar potential is given by 

V(h, H, A, if±) = ^ [ 2 n\{h + v ) 2 + Ai(/i + v ) 4 + 2 /x 2 (A 2 + i / 2 + 2 H+H~) 

+X 2 (A 2 + H 2 + 2H + H~) 2 ] 

+\{h + t -) 2 [A 3 H + H- + A S A 2 + A l H 2 ] 

where 

Al,s — 2 (A 3 + A 4 ± A 5 ) . 

Masses of these scalars are given by 

Ml = nl + 3Ain 2 , 

— ^2 + Az,n 2 , 

= ^2 + Asv 2 , 

Mjj± = nl + -A 3 n 2 . 

For A 4 — A 5 < 0 and A 5 > 0 (A 4 + A 5 < 0 and A 5 < 0), A (H) is the lightest Z 2 odd particle 
(LOP). In this work, we take A as the LOP and hence, as a viable DM candidate. Choice of H as 
LOP will lead to similar results. 

As we will explain later, for large DM mass, Ma 3> M z , appropriate relic- density of DM is 
obtained if Ma, Mh, and Mh± are nearly degenerate. Hence, in anticipation, we define 

AM h = Mh — Ma , 

AMh± = -D// : — Ma • 

so that the new independent parameters for the ID model become {Ma, AMh, AM h ±, A 2 , A 5 }. 
Here we have chosen A 5 as we treat A as the DM particle. 

The one-loop effective potential for h in the MS scheme and the Landau gauge is given by 

VfM +,D (ft) = n SM (ft) + Vj ID (h) (2.4) 


( 2 . 2 ) 

(2.3) 


where 


V^ih) = 


Hi 


5 

64tj-2 

2=1 


Mf{h) 


In 


MM. 

/i 2 (t) 


(2.5) 


n t is the number of degrees of freedom. For scalars and gauge bosons, n l is positive, whereas for 
fermions it is negative. Here Ch,G,f — 3/2, Cw.z = 5/6, and the running energy scale /j is expressed 
in terms of a dimensionless parameter t as /i(t) = M z exp(f). M, is given as 


Mi{h) — Ki(t) h 2 (t) — «/(£). 
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n l , Ki and k/ ; can be found in Eq. (4) in Ref. [21] (see also Refs. 

The additional contribution to the one-loop effective potential due to the inert doublet is given 

by [32! 


Vi lD (h) = 


1 


647T 2 

j=H,A,H+,H~ 


MHh) 


In 


Mf(h) 

/i 2 (f) 


( 2 . 6 ) 


where 

Mf(k)= 1 -X j (t)h 2 (t) + I 4(t) (2.7) 

with A A (t) = 2A s (t), A H (t) = 2A L (f), and A^±(f) = A 3 (f). 

In the present work, in the Higgs effective potential, SM contributions are taken at the two-loop 
level 0 0 12SI EZj, whereas the ID scalar contributions are considered at one loop only. 

For h v , the Higgs effective potential can be approximated as 


with 


where [6] 


Ai,,ff(fc) = 0*0 + A”,(A), 


( 2 . 8 ) 

(2.9) 


Af“«(A) = e 4r(W |A t (/j = h) + A^„(fi = h) + A^(/j = ft) 


,P> 


Aj D « a (/0= E 


=4T(fc) 


j=L,S,3 


64 ^ 1111( ' 5 ^ ) - 2 


Here 5,- = 1 when j — L, S] Sj — \ for j — 3; and 


f h 

T (h) — / y(/i) d ln/^. 

JM t 


Anomalous dimension y(/i) of the Higgs field takes care of its wave function renormalization. As 
quartic scalar interactions do not contribute to wave function renormalization at the one-loop level, 
ID does not alter 7 (/r) of the SM. The expressions for the one- and two-loop quantum corrections 
A) in the SM can be found in Ref. [6|. All running coupling constants are evaluated at [i = h. 
This choice guarantees that the perturbative expansion of the potential is more reliable [25], which 
can be understood as follows. The perturbative expansion of the potential at high field values of 
h consists of terms of the form e In (h/n), where e is some dimensionless coupling. Clearly the 
perturbation series works better if we choose the renormalization scale /j ~ h and ensure that e 
is small as well. Hence, our choice of n — h, along with the requirement that all the couplings 
remain within the perturbative domain, ensures faster convergence of the perturbation series of 
the potential. 
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To compute the RG evolution of all the couplings, we first calculate all the SM couplings at M t , 
taking care of the threshold corrections as in Refs. [41IT1 18l 129]. Then we evolve them up to M-p\ using 
our own computer codes incorporating the RG equations [301433] . The corresponding /3-functions 
include three-loop SM effects and two-loop ID contributions (one-loop RGEs are presented in the 
Appendix for completeness). If DM mass Ma is larger than M t , then the ID starts to contribute 
after the energy scale Ma- For Ma < M t , the contributions of ID to the /3-functions are rather 
negligible for the running from Ma to M t , as is evident from the expressions. 

To help the reader in reproducing our results, we provide in Table [T] a specific set of values of 
A i at M t = 173.1 GeV and at M P1 = 1.2 x 10 19 GeV for M h = 125.7 GeV and a s ( M z ) = 0.1184. 
In Fig. [I] we explicitly show running of the scalar couplings (Ai) for this set of parameters. We 
see that for this specific choice of parameters, Ai assumes a small negative valu^] leading to a 
metastable EW vacuum as discussed in the following sections. This set is chosen to reproduce the 
DM relic density in the right ballpark. 



A s 

A L 

A 2 

A 3 

a 4 

A 5 

Ai 

M t 

0.001 

0.039 

0.10 

0.0399 

0.00003 

0.038 

0.127 

Mpi 

0.046 

0.082 

0.127 

0.090 

0.038 

0.036 

-0.009 


TABLE I. A set of values of all ID model coupling constants at Mt and Mp\ for Ma = 573 GeV, AM#± = 1 
GeV, A Mh = 2 GeV, and A 5 {Mz) = 0.001. 


* Ai and Ai. e ff, being almost equal at Mpi, are interchangeable at that scale. 
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FIG. 1. RG evolution of the couplings A; (i = 1, ..,5), A l, A s for the set of parameters in Table [IJ 

III. CONSTRAINTS ON ID MODEL: A BRIEF REVIEW 

ID model parameter space is constrained from theoretical considerations like absolute vacuum 
stability, perturbativity and unitarity of the scattering matrix. Electroweak precision measure¬ 
ments and direct search limits at LEP put severe restrictions on the model. The recent measure¬ 
ments of Higgs decay width at the LHC put additional constraints. The requirement that the ID 
DM saturates the DM relic density all alone restricts the allowed parameter space considerably. 
Although these bounds are already discussed in the literature [31], here we apply these bounds 
with the requirement that the model needs to be valid till M-p\. 


A. Vacuum stability bounds 

The tree level scalar potential potential U(<l>i, $ 2 ) is stable and bounded from below if [5] 

•^ 1 , 2 (A) > 0, As(A) > — 2 a/Ai(A)A 2 (A), Az, : s(A) > —a/Ai(A)A 2 (A) (3-1) 

where the coupling constants are evaluated at a scale A using RG equations. However, these 
conditions become nonfunctional if Ai becomes negative at some energy scale to render the EW 
vacuum metastable. Under such circumstances we need to handle metastability constraints on the 
potential differently, which we pursue in the next section. 
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B. Perturbativity bounds 


The radiatively improved scalar potential remains perturbative by requiring that all quartic 
couplings of V ($ 1 , $ 2 ) satisfy 


Al,2,3,4,5 |< 4vr 


(3.2) 


C. Unitarity bounds 


Unitarity bounds on A* are obtained considering scalar-scalar, gauge boson-gauge boson, and 
scalar-gauge boson scatterings [35]. The constraints come from the eigenvalues of the corresponding 
S-matrix [55] : 


A 3 ± A 

OO 

VI 

a 3 ± a 5 

< 

8:/r 


A3 - 

P 2A4 ± 3A5 

< 

8:/r 

— Ai — A2 

± \J (Ai 

— A2) 2 + A 2 

< 

87 T 

3 A 2 ± v / 9 (AT 

- A 2 ) 2 + 

(2A3 + A4) 2 

< 

8:/r 

— Ai — A2 

± yj (Ai 

- A 2 ) 2 + A 2 

< 

8:/r 


(3.3) 


D. Bounds from electroweak precision experiments 


Bounds ensuing from electroweak precision experiments are imposed on new physics models via 
Peskin-Takeuchi [37] S, T , U parameters. The additional contributions from ID are given by 


AS = — 

2tt 


1 , (Ml 

- In ' H 


5 

- — + 


MjjMl 


+ 


M\{M\ - 3 M- 


6 \M 2 h± J 36 3 (Ml-Mjj) 2 6 (M| - M 2 ) 3 


Hl In ( M 

Ml 


and 


AT = 


1 


327T 2 av 2 

where the function F is given by 


F + F ( Ml ± ,Ml ) - F (Ml M%) 




F(x,y ) — ( 2 x y 

0, x = y 

We use the NNLO global electroweak fit results obtained by the Gfitter group 

AS = 0.06 ±0.09, AT = 0.1 ±0.07 


(3.4) 


(3.5) 


(3.6) 


(3.7) 














> 

o 

% 

<1 



50 100 150 200 

AM# [GeV] 


FIG. 2. Allowed parameter space in A M H ± — A Mu plane for Ma = 70 GeV and A g = 0.007. Constraints 
from S and T parameters are shown by solid black, green and red lines. The blue band corresponds the 
3a variation in f lh 2 = 0.1198 ± 0.0026 jT0]. On the brown region the unitarity bound is violated on or 
before Mp\. The cross-hatched region is excluded from LEPII data. 


with a correlation coefficient of +0.91, fixing AU to zero. We use this fit result as the contribution 
of the scalars in the ID model to AU is rather negligible. 


To assess the implications of S and T constraints on the ID model, in Fig. [2j in the AM H ± —AM H 
plane, we display various constraints for Ma = 70 GeV and A g = 0.007. This is the maximum 
value of Xs for the given DM mass, allowed by LUX HU direct detection data at la. The blue 
region allowed by relic density constraints shifts upwards for smaller A g. Constraints on AS and 
AT, as mentioned in Eq. (3.7), are marked as black, green, and red solid lines. We see that the la 
bound on AT is the most stringent one. The black line corresponding to the lower limit on AS at 
la can also be interesting. But the LEP II bounds, represented by the cross-hatched bar, take away 
a considerable part. The line representing the AS upper limit is beyond the region considered and 
lies towards the bottom-right corner of the plot. In this plot, increasing Mh enhances S, but T 
gets reduced. In addition, if we demand unitarity constraints to be respected up to Mpi, assuming 
no other new physics shows up in between, the parameter space gets severely restricted. I 11 this 
plot, a small window is allowed by AT only at 2a, which satisfies DM relic density constraints. As 
mentioned earlier, this window gets further reduced with smaller A# (the relic density band takes 
an “L” shape as shown in Fig. [5]). 
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E. Direct search limits from LEP 


The decays Z —* AH, Z —> H + H~, W ± —>• AH ± , and W ± —> HH ± are restricted from Z and 
W ± decay widths at LEP. It implies Ma + Mh > Mz, 2 M H ± > Mz , and M H ± + Mh,a > 

More constraints on the ID model can be extracted from chargino [42! and neutralino [53] search 
results at LEP II: The charged Higgs mass M H ± > 70 GeV. The bound on Ma is rather involved: 
If Ma < 80 GeV, then A M H should be less than ~ 8 GeV, or else M H should be greater than 
~ 110 GeV (see Fig. [4]). 


F. Bounds from LHC 


In the ID model, Higgs to diphoton signal strength is defined as 

_ a(gg ->•/&->• 77 ) _ Br(h -)• 77 ) ID 

^ 77 a(gg h 77 )sm Br(h 77 )sm 


(3.8) 


using the narrow width approximation for the production cross section of cr(gg —)■ h —)■ 77 ) and 
the fact that <j(gg —> h ) in both the SM and ID are the same. 

Now if the ID particles have masses less than Mh/ 2, h —y ID, ID decays are allowed. In that 
case, 


r(/i —>■ 77)10 r tot (/i^SM,SM) 

r(ft -> 77 )s M r M (ft -> sm, sm) + r to ,(h -> id, id) ’ 


(3.9) 


where [55] 


r (h-> ID, ID) 


V 2 ,2 7 4Mg,V /2 

1 6nM h 10 V Ml ) 


(3.10) 


where for ID = A, H, H ± , Aid = As, Al, y/2 A 3 . 

In this case, the ID particles are heavier than Mh/ 2, 


_ n h -»• 77)id 
f(ft -> 77) S M 


(3.11) 


In the ID model, the H ± gives additional contributions at one loop. The analytical expression is 
given by ]5ol [57] 


r (h -)■ 77 ) id = 


am 


h 


2567r 3 r> 2 


^ NfQ 2 fyfFi/ 2 (rf) + ywF\{r w ) + Qh± 


Vy-hH+H- 
2 m 2 H± 


Fq(t h ±) 


(3.12) 


where r t = m\/Arrij. Qf, Qh± denote electric charges of corresponding particles. N/ is the color 
factor, i/f and yw denote Higgs couplings to ff and WW. yhH+H- — A 3 V stands for the coupling 
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constant of the hH + H vertex. The loop functions ^(o, 1 / 2 , 1 ) are defined as 

F o{r) = -[t-/(t)]t~ 2 , 

F 1/2 (t) = 2[r + (r - l)/(r)]W 2 , 

Fi(r) = — [2r 2 + 3r + 3(2r - l)/(r)]W 2 , 


where 


f(r) = 


Sill 


-1 




■sK 


1 | y 7 1 — r 1 
\/l — T~ 1 


Z7T 


r < 1 
r > 1 


(3.13) 


From the diphoton decay channel of the Higgs at the LHC, the measured values are /r 77 =1.17 ± 
0.27 from ATLAS [3BJ and /i 77 = 1 . 141 q ; 2 3 from CMS [39] . 

One can see that a positive A 3 leads to a destructive interference between SM and ID contribu¬ 
tions in En. (3.12) and vice versa. Hence, for ID particles heavier than Mh/ 2, /x 77 < 1 (/i 77 > 1) 
when A 3 is positive (negative). However, if these ID particles happen to be lighter than M h /2, 
they might contribute to the invisible decay of the Higgs boson. Using the global fit result [SOJ 
that such an invisible branching ratio is less than ~ 20%, in Eq. (3.9), the second ratio provides a 
suppression of ~ 0.8 - 1 . 

Now can we work with a negative A 3 in the ID model? We will discuss this at the end of this 
section. For the benchmark points used in this paper, we have worked only with positive values 
of A 3 , allowed at lcr by both CMS and ATLAS experiments. It is important to note that the 
constraints from /r 77 pass through the sign of A 3 in our vacuum stability considerations. 


G. Constraints from dark matter relic density and direct search limits 

The ID dark matter candidate A can self-annihilate into SM fermions. Once the DM mass is 
greater than the IE-mass, so that the DM can annihilate into a pair of W bosons, the cross section 
increases significantly, thereby reducing DM relic density. Hence, it becomes difficult to saturate 
fill 2 after rs./ 75 GeV with positive A 5, although for M4 < 75 GeV, both signs of A s can be allowed 
to arrive at the right DM relic density flh 2 . 

The role of the sign of A 5 can be understood from the contributing diagrams to the AA —y 
W + W~ annihilation processes. Four diagrams contribute: the AAW + W~ vertex driven point 
interaction diagram (henceforth referred to as the p-channel diagram), /7 + -mediated t- and u- 
channel diagrams, and the h- mediated s-channel diagram. For A A —yZZ annihilation, the t- and 
n-channel diagrams are mediated by H. A negative A 5 induces a destructive interference between 
the s-channel diagram with the rest, thereby suppressing AA —> W + W~, ZZ processes. For DM 
masses of 75 — 100 GeV, this can be used to get the appropriate Q/ 1 2 [ 12 ], 1ST] . To avoid large 
contributions from t- and w-channel diagrams and coannihilation diagrams, the Mh and M#± can 
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be pushed to be rather large > 500 GeV. However, to partially compensate the remaining p-channel 
diagram by the s-channel one, A 5 assumes a large negative value ~ — 0 . 1 , which is ruled out by 
the DM direct detection experiments. That is why in the ID model, DM can be realized below 75 
GeV, a regime we designate as the “low” DM mass region. 

At “high” DM mass Ma > 500 GeV, one can get the right Vth 2 due to a partial cancellation be¬ 
tween different diagrams contributing to the AA —>■ W + W~ and AA —>■ ZZ annihilation processes. 
For example, in AA —)■ W + W~, the p-charmel diagram tends to cancel with the // + -mediated /- 
and u-channel diagrams m E 2 ] in the limit Ma Mw, and the sum of amplitudes of these 
diagrams in this limit is proportional to M^± — M Hence, at high Ma, a partial cancellation 
between these diagrams is expected for nearly degenerate M H ± and Ma- Similarly, for AA —s- ZZ, 
a cancellation is possible when the masses Mh and Ma are close by. For Ma > 500 GeV, keeping 
the mass differences of M#± and Mh with Ma within 8 GeV, such cancellations help reproduce 
the correct Vth 2 . It is nevertheless worth mentioning that such nearly degenerate masses will lead 
to coannihilation of these Z 2 odd ID scalars [53] to SM particles. Despite such near degeneracy, 
both H and H + , being charged under the same Z 2 as A, decay promptly to the LOP A, so that 
they do not become relics. We use FeynRules [53] along with micrOMEGAs [53156] to compute the 
relic density of A. 

DM direct detection experiments involve the /i-mediated t-channel process AN —>■ AN with a 
cross section proportional to \ 2 S /M\ in the limit Ma Mat: 


&a,n 


2 

m r 

71 



Ag 

M a M 2 


2 


(3.14) 


where / ~ 0.3 is the form factor of the nucleus. m r represents the reduced mass of the nucleus 
and the scattered dark matter particle. 

Thus, A 5 is constrained from nonobservation of DM signals at XENON 100 [57, 58] and LUX BU- 
For Ma = 70 GeV, the ensuing bound from LUX jHj, [59] data at lcr is |Ag| < 0.007. 

The constraint on A 5 from DM direct detection experiments gets diluted with Ma [see 
Eq. (3.14)]. Hence, for low DM mass, direct detection bounds are more effective. At high mass, the 
relic density constraints are likely to supersede these bounds. For example, for Mom = 573 GeV, 
the upper limit on |Ag| is 0.138 from LUX. However, to satisfy the relic density constraints from 
the combined data of WMAP and Planck within 3<r, A s can be as large as 0.07 only. 

Within the framework of the ID model, it is possible to explain the observations in various 
indirect DM detection experiments [15, !60] for some regions of the parameter space. In this paper, 
however, we do not delve into such details as such estimations involve proper understanding of 
the astrophysical backgrounds and an assumption of the DM halo profile which contain some 
arbitrariness. For a review of constraints on the ID model from astrophysical considerations see, 
for example, Ref. [61]. 
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Sign of A3 


Whether A 3 can be taken as positive or negative depends on the following: 

• If the ID model is not the answer to the DM puzzle, so that both relic density and direct 
detection constraints can be evaded, no restriction exists on the possible sign of A 3 (M Z ). 
Otherwise, the following two cases need be considered: 


— A negative X 3 (M Z ) implies 

A s {M z )<--^M 2 h± -M 2 a ). 
v 

As we are considering A as the DM candidate, so that M a < M H ±, A s(M z ) is always 
negative when X 3 {M Z ) < 0 . For low DM mass, the splitting ( M H ± — M A ) > 10 GeV, 
as otherwise DM coannihilation processes cause an inappreciable depletion in Q/ 1 2 . For 
M a = 70 GeV, this implies a lower bound A s(M z ) < —0.025, which violates the DM 
direct detection bound |As[ < 0.007. Hence, for low DM mass, a negative A 3 (M Z ) is 
not feasible. 

— For high mass DM, the right relic density can be obtained when the splitting ( Mh± — 
M a ) ~ a few GeV or less. The above logic then implies that a negative A 3 does not put 
any severe restriction on A 5 to contradict DM direct detection bounds as earlier. Hence, 
for a high DM mass, A 3 {M Z ) can assume both the signs. Moreover, due to propagator 
suppression for large M#± in the /177 vertex, the ID contribution to /x 77 is negligibly 
small and hence, the sign of A 3 (M Z ) is not constrained by measurements on yU 77 as well. 


If at any scale, A 3 is negative while Ai > 0, then the bound (3.1) must be respected. 


• If at some scale, Ai < 0, then a negative A 3 makes the potential unbounded from below, as 
mentioned in the following section. This means one can start with a negative X 3 (M Z ), but 
with RG evolution when Ai turns negative, A 3 evaluated at that scale must be positive. Such 
parameter space does exist. Here, we note a significant deviation of our analysis from earlier 
analyses which did not allow a negative Ai. For example, in Ref. [0] a negative A 3 (M z ) was 
not allowed from stability of the Higgs potential if the theory has to be valid up to 10 16 GeV 
together with relic density considerations. 


IV. TUNNELING PROBABILITY AND METASTABILITY 

In the standard model, the present measurements on M\ and M t indicate that the electroweak 
vacuum, in which the Universe is at present residing, may be a false one. From this metastable 
vacuum, it might tunnel into a deeper true vacuum, residing close to Mp\. 
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(a) 


(b) 


FIG. 3. (a) Tunneling probability Vo dependence on M t . The left band (between dashed lines) 

corresponds to SM. The right one (between dotted lines) is for ID model for DM mass Ma = 573 GeV. 
Constraints from WMAP and Planck measured relic density, as well as XENON 100 and LUX DM direct 
detection null results are respected for these specific choice of parameters. Light-green band stands for 
Mt at ±l(j. (b) Vo is plotted against Higgs dark matter coupling A s(Mz) for different values of A 2 (Mz)- 


The decay probability of the EW vacuum to the true vacuum at the present epoch can be 
expressed as [6l [62], [63] 

A 4 

(4.1) 


Vo = 0.15|| e- s(Afl) 
where H is the Hubble constant and the action is given by 

„,. , 8tt 2 

S(A b ) = 


3|Ai(A s )| 


(4.2) 


S(Ab) is called the action of bounce of size R = A^ . The value of R for which S(Ab) is min¬ 
imum gives the dominant contribution to the tunneling probability Vo■ It occurs when Ai(A^) 
is minimum, i.e., /3\ 1 (Ab) = 0. Henceforth, A b denotes the scale where Ai is minimum. Here 
we neglect the loop correction to the action, as in Ref. [62] it had been argued that setting the 
running scale to R significantly restricts the size of such corrections. In the numerical evaluation 
we have used Ai je ff in place of Ai in Eq. (4.2). We also neglect gravitational corrections [Ml 165] to 
the action as in Ref. [4]. In Ref. |62j it was pointed out that thermal corrections are important 
at very high temperatures. Finite temperature effects to EW vacuum stability in the context of 
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SM have been calculated recently in a preprint It had been claimed that the parameter 

space corresponding to EW metastability shrinks considerably. In this paper we work with zero 
temperature field theory only. 

The presence of the inert doublet induces additional contributions to f3\ l [see Eq. (A2)]. As 
a result, which is generic for all scalars, Ai receives a positive contribution compared to the SM, 
which pushes a metastable vacuum towards stability, implying a lower Vo- 

Electroweak metastability in the ID model has been explored earlier in the literature, albeit in a 
different context [68HTT] . If H + gets a VEV, there could exist another charge-violating minimum. 
If instead, A receives a VEV, another (7P-violating minimum could pop up. But these vacua 
always lie higher than the usual EW vacuum. If Z 2 is broken by introducing additional soft terms 
in the Lagrangian, then the new Z 2 -violating minimum can be lower than the usual Z 2 -preserving 
EW minimum. As in our present work, Z 2 is an exact symmetry of the scalar potential, so such 
cases need not be considered. However, as mentioned earlier, if at some scale before Mpi, the sign of 
Ai becomes negative, there might exist a deeper minimum which is charge-, CP- and Z 2 -preserving 
and lying in the SM Higgs h direction. 

Whether the EW vacuum is metastable or unstable, depends on the minimum value of Ai before 
Mpi, which can be understood as follows. For EW vacuum metastability, the decay lifetime should 
be greater than the lifetime of the Universe, implying Vo < 1. This implies [4, (721 

-0.06488 


Al,eff(A_B) > Ai ,min (A b ) — 


1 - 0.00986 In {v/K b) 


(4.3) 


Hence we can now reframe the vacuum stability constraints on the ID model, when Ai ie ff runs into 
negative values, implying metastability of the EW vacuum. We remind the reader that in the ID 
model, instability of the EW vacuum cannot be realized as addition of the scalars only improves 
the stability of the vacuum. 


• If 0 > Ai ie fr(As) > Ai m i n (Ap}), then the vacuum is metastable. 

• If Ai i eff(A j B) < Ai )min (A B ), then the vacuum is unstable. 

• If A 2 < 0, then the potential is unbounded from below along the H , A and H ± direction. 

• If As(Ai) < 0, the potential is unbounded from below along a direction in between H ± and 
h. 

• If A jL (A i ) < 0, the potential is unbounded from below along a direction in between H and h. 

• If A s (A 0 < 0, the potential is unbounded from below along a direction in between A and h. 


I See Ref. [67] for an earlier work. 
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In the above, A T represents any energy scale for which Ai )e fj is negative and the conditions for 
unboundedness of the potential follow from Eq. ( 2 . 2 ). At this point note the significant deviations 
we are making in the allowed parameter space compared to the usual vacuum stability conditions: 
According to Eq. (3.1), A 3^,5 can take slightly negative values. But at a scale where Ai >e ff is 
negative, with the new conditions A 3i l,s have to be positive. 


The tunneling probability Vq is computed by putting the minimum value of Ai e ff in Eq. (4.2) to 
minimize S'(As). In Fig. |3](a), we have plotted Vo as a function of M t . The right band corresponds 
to the tunneling probability for our benchmark point as in Table [TJ For comparison, we plot Vo 
for SM as the left band in Fig. |3](a). la error bands in a s and M^ are also shown. The error due 
to a s is more significant than the same due to M^. As expected, for a given M t , the presence of 
ID lowers tunneling probability. This is also reflected in Fig. [3j(b) , where we plot Vo as a function 
of A s(Mz) for different choices of A 2 (Mz), assuming Mh = 125.7 GeV, M t = 173.1 GeV, and 
a s = 0.1184. Here DM mass M A is also varied with X s to get Xlh 2 = 0.1198. For a given X S (M Z ), 
the higher the value of A 2 (M Z ), the smaller Vq gets, leading to a more stable EW vacuum. 


V. PHASE DIAGRAMS 

The stability of EW vacuum depends on the value of parameters at low scale, chosen to be M z . 
In order to show the explicit dependence of EW stability on various parameters, it is customary 
to present phase diagrams in various parameter spaces. 

In Fig. [4] we show the LEP constraints in the Mh — M A plane as in Ref. [43]. We update this 
plot identifying regions of EW stability and metastability. As we are considering a scenario where 
the ID model is valid till Mp\, there are further limits from unitarity. The relic density constraint 
imposed by WMAP and Planck combined data is represented by the thin blue band. The choice 
of A 2 (M z ) does not have any impact on relic density calculations, but affects EW stability as 
expected. In this plot, for higher values of A 2 (Mz), the region corresponding to EW metastability 
will be smaller. The chosen parameters satisfy the LUX direct detection bound. 

As small splitting among M A , Mh, and Mh± leads to some cancellations among diagrams 
contributing to DM annihilation, A M#± and A Mh are often used as free parameters in the ID 
model. In Fig. [5j we present constraints on this parameter space for M A = 70 GeV. As before, the 
brown region corresponds to unitarity violation before Mp\. For small A M H ± and A Mh, As^s are 
required to be small, which leads to little deviation from SM metastability. The metastable region 
is shown by the yellow patch, which shrinks for larger A 2 . The blue band reflects the relic density 
constraint for A s{M z ) = 0.001. For such small A s{M z ), the h-mediated s-channel diagram in 
AA —y WW or AA —> ZZ contributes very little. H + - or //-mediated t- and w-channel diagrams 
are also less important than the quartic vertex driven diagram due to propagator suppression. 
This explains the “L” shape of the blue band. For higher values of A s{M z ), the shape of the band 
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FIG. 4. Constraints in Mh — Ma plane. The cross-hatched region is excluded from LEP jl3j. Choosing 
M h ± = 120 GeV and A s{Mz) = 0.001, relic density constraint is satisfied at 3c on the blue band. 
The green (yellow) region corresponds to EW vacuum stability (metastability). The solid brown line 
correspond to Mh = Ma ■ The grey area on the left to it is of no interest to us as we have chosen 
Mh > Ma ■ The dashed brown line shows the LEP I limit. On the brown region, unitarity constraints are 
violated before Mp\. 


changes and ultimately leads to a closed contour. It appears that due to LEP constraints, EW 
vacuum metastability is almost ruled out. Although the LEP constraint permits A Mh < 8 GeV, 
allowing a narrow strip towards the left, the relic density constraints cannot be satisfied on this 
strip as it leads to an increased rate of DM coannihilation processes, leading to a dip in Qh 2 . But 
as we will see later, if M t and a s are allowed to deviate from their respective central values, for 
some region in this parameter space, it is possible to realize a metastable EW vacuum. 

To delineate the role of Mh in EW vacuum stability, in Figs. [4] and [5j A 2 (Mz) was chosen to 
be small. Now to demonstrate the effect of As, we will now present in Fig. [6] phase diagrams 
in the A s{Mz) — Ma plane. Panel (a) deals with low DM masses. For A M H ± = 40 GeV and 
A Mh = 40 GeV, part of the allowed relic density band (blue) is allowed from LEP constraints 
(cross-hatched band). The entire parameter space corresponds to EW vacuum stability. Choosing 
small AM H ± and A Mh, which imply small values of As^s, can lead to metastability. But those 
regions are excluded by LEP. Again, metastability can creep in if M t and a s are allowed to deviate 
from their central values. 

In Fig. [6](b), we study the same parameter space for high DM masses. As mentioned before, 
to obtain the correct relic density, smaller mass splitting among various ID scalars needs to be 
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FIG. 5. Phase diagram in A Mh — A Mjj± plane for Ma = 70 GeV. The green and yellow regions 
correspond to EW vacuum stability and metastability respectively. The cross-hatched band is excluded 
from LEP. The brown region suffers from unitarity violation before Mp\. The blue band reflects relic 
density constraint at 3a. 


chosen. For AM#± = 1 GeV and A Mh = 2 GeV, the 3cr relic density constraint is shown as the 
bine band. The blue dashed line demarcates the boundary between stable (green) and metastable 
(yellow) phases of EW vacuum. The choice of small values of A M H ± and A M H , in turn, leads to 
a large region pertaining to EW metastability. 

To illustrate the sensitivity to the mass splitting, in Fig. |6](b), we present another relic density 
band (red) when A Mh± = 5 GeV and A Mh = 2 GeV. The corresponding boundary between the 
phases is denoted by the red dashed line. The region on the right implies EW stability (the green 
and yellow regions do not apply to this case). As for high DM masses, EW metastability can be 
attained for a sizable amount of the parameter space; A 2 (Mz) need not be chosen to be very small 
to maximize the metastable region for the sake of demonstration. 

The fact that for SM the EW vacuum stability is ruled out at ~ 3a, is demonstrated by a 
phase space diagram in the M t — M h plane [HI d- In Ref. [3], similar diagrams were presented 
for a singlet scalar extended SM. To demonstrate the impact of ID scalars to uplift the EW 
vacuum metastability, we present phase diagrams in the M t — Mh plane for two sets of benchmark 
points in Fig. [7J Panel (a) is drawn for Ma = 70 GeV, A M H ± = 11.8 GeV, A M H = 45 GeV, 
A s{M z ) = 0.001, and A 2 (Mz) = 0.1. For panel (b) the set of parameters in Table [T] is being used. 
Both sets of parameters are chosen so that they respect the WMAP and Planck combined results 
on DM relic density and the direct detection bounds from XENON 100 and LUX. As in Ref. [3], the 
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A s(M z ) 


(a) 


(b) 


FIG. 6. Phase diagram in A s(Mz) — Ma plane for A 2 (M z ) = 0.1. Panel (a) stands for ‘low’ DM mass. 
The blue band corresponds the 3cr variation in Qh 2 when A M H ± = 40 GeV and A Mu = 40 GeV. LEP 
direct search constraints are represented by the cross-hatched band at the bottom. Entire green region 
imply EW vacuum stability. Panel (b) stands for ‘high’ DM masses. The relic density band (blue) now 
correspond to A M H ± = 1 GeV and A Mh = 2 GeV. The corresponding stable and metastable phases 
for EW vacuum are represented by green and yellow patches respectively. The relic density band (red) 
corresponds to A M H ± = 5 GeV and A Mh = 2 GeV. For this, the boundary separating the EW phases 
is denoted by the red dashed line. 


line demarcating the boundary between stable and metastable phases of EW vacuum is obtained 
by demanding that the two vacua be at the same depth, implying Ai(As) = /3 a, (Ar) = 0. The line 
separating the metastable phase from the unstable one is drawn using the conditions /3\ 1 (A b ) = 0 
and Ai(A s ) = Ai ^J Ar), as in Eq. (4.3). The variations due to uncertainty in the measurement 
of a s are marked as dotted red lines. In each panel, the dot representing central values for Mh and 
M t is encircled by la, 2a, and 3a ellipses representing errors in their measurements. According 
to Fig. [7](a), EW vacuum stability is allowed at 1.5cr, whereas in Fig. [7](b), it is excluded at 2.1 a, 
indicated by blue-dashed ellipses. 

As in the literature SM EW phase diagrams are also presented in the a s (M z ) — M t plane [8, 72) . 
we do the same in the ID model as well. I 11 Fig. [8j we use the same sets of benchmark parameters 
as in Fig. [7J As a consistency check, one can note that the EW vacuum is allowed or ruled out at 
the same confidence levels. 

To study the impact of nonzero ID couplings, however, it is instructive to study the change in 
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A S (M Z )=0.001, M dm =70 GeV A S (M Z )=0.001, M dm =573 GeV 




(a) 


(b) 


FIG. 7. Phase diagrams in Mh — Mt plane. Panels (a) and (b) stand for ‘low’ and ‘high’ DM masses 
respectively. Regions of absolute stability (green), metastability (yellow), instability (red) of the EW 
vacuum are also marked. The grey zones represent error ellipses at 1, 2 and 3 a. The three boundary lines 
(dotted, solid and dotted red) correspond to a s [Mz ) = 0.1184 ± 0.0007. Details of benchmark points are 
available in the text. 


the confidence level (cr) at which EW stability is modified with respect to these couplings. As in 
Ref. pE], we plot in Fig. [9] a against Xs(Mz) for different values of A 2 (Mz). We vary Ma along with 
Xs(M z ) to keep DM relic density fixed at Qh 2 = 0.1198 throughout the plot. Note that changing 
X 2 (M z ) does not alter Vth 2 . The masses of other ID particles are determined using A Mjj± = 1 GeV 
and AMh = 2 GeV. The parameter space considered does not yield too large DM-nucleon cross 
section, inconsistent with XENON 100 and LUX DM direct detection null results. For a specific 
value of X 2 (Mz) = 0.1, with the increase of Xs(Mz), the confidence level at which EW is metastable 
(yellow region) gets reduced and becomes zero at Xs(M z ) — 0.04. After this, EW vacuum enters 
in the stable phase (green). With further increases in As(Mz), the confidence level at which EW 
is stable keeps increasing. To illustrate the role of X 2 (Mz), we use two other values in the same 
plot. The value of Xg(Mz) at which the EW vacuum enters in the stable phase increases with 
decreases in X 2 (Mz), as expected. The yellow and green marked regions are not applicable when 
X 2 [M z ) = 0.05,0.15. 
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A 5 (M z )=0.001, M dm =70 GeV 



(a) 


FIG. 8. Phase diagrams in Mt — a s (Mz ) plane 
Notations used are also the same as in Fig. [7J 



172 174 176 178 180 


Top pole mass M, [GeV] 

(b) 

the same sets of benchmark points as in Fig. [7| 



FIG. 9. Dependence of confidence level at which EW vacuum stability is excluded (one-sided) or allowed 
on A s{Mz) and A 2 {Mz)- Regions of absolute stability (green) and metastability (yellow) of EW vacuum 
are shown for A 2 (Mz) = 0.1. The positive slope of the line corresponds to the stable electroweak vacuum 
and negative slope corresponds to the metastability. 
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VI. VELTMAN’S CONDITIONS 


As we have extended the validity of the ID model till Mp\, it is interesting to explore whether 
Veltman’s condition (VC) can be satisfied in this model at any scale on or before Mp\. It is 
particularly interesting as similar studies for SM have been carried out in Ref. [7]. It has been shown 
that if one imposes VC in the SM at Mp\, then the top mass measurement M t = 173.1 ± 0.6 GeV 
implies Mh ~ 135 ± 2.5 GeV, which is excluded at more than 3cr. 

Veltman’s condition implies that the quadratic divergences in the radiative corrections to the 
Higgs mass can be handled if the coefficient multiplying the divergence somehow vanishes [731 El]. 
VC includes the contributions from the infrared degrees of freedom of the theory and does not carry 
any special information about the ultraviolet divergences. In SM, it suggests the combination 

6 Ai + -gl + -<7i — 12 y 2 t = 0. 

Due to the large negative contribution from the term containing the top Yukawa coupling, it is 
not possible to satisfy VC till Mp\ given the experimental measurements of M t and M\ within the 
context of SM. 

In the ID model, as we add more scalars a possibility opens up to satisfy VC, as their contri¬ 
butions can offset the large negative contribution from the top quark. 

The above VC for the SM associated with /ii is promoted in the ID model to [38, 75] 

6Ai + 2A3 + A4 + -g 2 + -gl — 12 yl = 0. (6.1) 

If 2 A 3 + A 4 is positive, then we have checked with our RG improved coupling constants, and it is 
possible to satisfy the above VC at a scale before M P1 . 

However, in the ID model /i 2 also receives quadratically divergent radiative corrections. The 
corresponding VC reads as 

6 A 2 + 2 A 3 + A 4 + -gl + -gf = 0 . ( 6 . 2 ) 

Note that it lacks the Yukawa contribution as the unbroken Z 2 forbids fermionic interactions of 
the inert doublet $ 2 . As 2A 3 + A 4 is already positive, this VC can be satisfied if A 2 is negative. But 
a negative A 2 renders the potential unbounded from below as evident from our earlier discussions. 
Note that amongst our RG improved coupling constants Ai can be driven to negative values at high 
scales. But this makes the required cancellations for VCs even worse. Hence, it is not possible to 
satisfy Veltman’s conditions in a scenario where only the ID model reigns the entire energy regime 
up to the Mp\. 
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VII. SUMMARY AND CONCLUSION 


If the standard model is valid up to the Planck scale, the present measurements on the masses 
of the top quark and Higgs indicate the presence of a deeper minimum of the scalar potential at 
a very high energy scale, threatening the stability of the present electroweak vacuum. State of 
the art NNLO calculations performed to evaluate the probability that the present EW vacuum 
will tunnel into the deeper vacuum lying close to Mp\ suggest that the present EW vacuum is 
metastable at ~ 3c. The lack of stability might be the artifact of incompleteness of the SM. 

Although the LHC has yet to find any signal suggesting existence of any new physics beyond 
the standard model of particle physics, other experimental evidence points towards the existence 
of dark matter, which so far could have escaped detection in colliders and DM direct detection 
experiments. Hence, it is important to look into the problem of EW stability in a scenario which 
addresses the issue of DM as well. In particular, we extend SM by adding an inert scalar doublet, 
offering a viable DM candidate and assume that this model is valid up to Mp\. 

In this paper, our intentions are twofold. First, in such a scenario we have consolidated the 
bounds imposed on the ID model. As we are demanding validity of the model up to Mpi, the RG 
evolution of the couplings can disturb the unitarity of the S-matrix governing various scattering 
processes, which in turn imposes stringent limits on the parameter space at the EW scale. In this 
light, we present a consolidated discussion updating the existing bounds on the ID model. 

The other and the main goal of this paper is to check the stability of EW vacuum in the ID model. 
If the ID DM happens to be the only DM particle, which saturates the observed DM relic density, 
can the ID model modify the stability of the EW vacuum? Note that rather than considering new 
physics effects close to the Planck scale, here new physics is added at the EW scale only. It is well 
known that addition of a scalar can improve stability of the EW vacuum. But if we are to solve 
both the DM and EW vacuum stability problems in the context of the ID model, it is important to 
study the parameter space which allows us to do so. As ID introduces a few new parameters and 
fields, the study of the parameter space is quite involved when we consider radiatively improved 
scalar potentials containing SM NNLO corrections and two-loop ID contributions. Inclusion of 
these NNLO corrections is mandatory to reproduce the correct confidence level at which EW 
vacuum is metastable in the SM. In case of ID, if one works with one-loop ID contributions in 
the RGEs instead of two-loop ones, the changes in the plots presented in the paper are limited 
to 5% only. However, as a small change in the action S is amplified exponentially in the vacuum 
decay lifetime, two-loop effects have also been taken care of. For the benchmark point used in 
Table [IJ Vq using one-loop ID RGEs turns out to be 5.8 x 10 -1271 , whereas for two-loop ID RGEs, 
Vo changes to 1.4 x 1CT 1310 . We neglect other loop corrections to bounce action, as a 5%-10% 
correction can lead to an appreciable change in the lifetime, but as demonstrated above with the 
example of inclusion of two-loop ID corrections, the phase diagrams essentially remain unchanged 


23 


(see Ref. [7B] for a related example). We note in passing that unlike Ref. p3j, our analysis is 
valid up to M P1 . Requiring the potential bounded from below at all scales below M P[ puts severe 
constraints on the model. 


As far as the allowed parameter space is concerned, we focus on metastability of EW vacuum 
which is realized when Ai is negative. In this context, our analysis is much different from earlier 
analyses in the literature, for example, as in Ref. m- When Ai is negative, the vacuum stability 
conditions (3.1) are no longer useful and the parameter space analysis in Ref. ra is not valid. 
For example, when Ai is positive, A 3^5 can assume small negative values as positivity of the Ai h 4 
term ensures the potential bounded from below. However as we point out, one can no longer 
afford a negative A 3^,5 at a scale where Ai is negative. If we demand in addition that the ID 
model saturates DM relic density, then even when Ai is positive, A 3 cannot be negative. But as 
mentioned earlier, for negative Ai, although A 3 has to be positive at that scale, its value at Mz can 
still be negative. Thus, one of the main motives for this work is to find out the allowed parameter 
space in the ID DM model in case the EW vacuum is metastable. 


We see that for DM masses of 70 GeV, the allowed parameter space corresponds to absolute 
stability unless we allow some deviation of M t and a s from their measured central values. For 
higher DM masses more than 500 GeV, it is possible to realize a metastable EW vacuum for a 
large parameter space. This vacuum will have a longer lifespan than the SM one as the addition 
of scalars improves the stability of the EW vacuum. 
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Appendix A: One-loop beta functions 


The beta functions of the quartic coupling parameters for the ID model are defined as 


P\i 


16tt 2 


d\j 
d In (i 


(AI) 
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The expressions at one loop are given by mm 

/3 a x — 24A| + 2A 1 2 3 4 5 6 7 8 + 2 A 3 A 4 + A 2 + A 2 

+ g (3^2 + 9i + ~ 3Ai (3y 2 + 9i) 

+4Ai (y 2 + 3y 2 + 3y 2 ) — 2 (y* + 3yj) + 3yf ) , (A 2 ) 

/3a 2 = 24A^ + 2Ag + 2 A 3 A 4 + A 4 + Ag 

+ g (3.92 + 9i + 2y 2 9i) — 3A 2 (3y 2 + y 2 ) , (A3) 

/3a 3 = 4 (Ai + A 2 ) ( 3 A 3 + A 4 ) + 4A 3 + 2A 2 + 2A 2 
+ ^ (392 + 9i ~ 2y29i) — 3A 3 (3y 2 + 9\) 

+2A 3 (y 2 + 3y 2 + 3y 2 ) , (A4) 

/3a 4 = 4 A 4 (Ai + A 2 + 2 A 3 + A 4 ) + 8 A 2 
+3y 2 9i — 3 A 4 (3y 2 + 9\) 

+ 2 A 4 (y 2 + 3y 2 + 3y 2 ) , (A5) 

/3a 5 = 4 A 5 (Ai + A 2 + 2 A 3 + 3 A 4 ) 

—3A 5 (3y 2 + y 2 ) 

+2A 5 (y 2 + 3y 2 + 3y 2 ) . (A 6 ) 

Let us note the y t dependence of these expressions. While /3\ 1 is dominated by the yf term, f3\ 2 
does not depend on y t . The y t dependence of other f3 \ i s are softened by the corresponding A* 
multiplying the y^ terms. Two-loop RGEs used in this work have been generated using SARAH [78j. 
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